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Abstract. We report the association of the INTEGRAL source IGR J 18027-20 16 with the BeppoSAX source SAX J 1802.7- 
2017. IGR J18027-2016 is seen to be a weak, persistent source by the IBIS/ISGRI instrument on board INTEGRAL with an 
average source count rate of 0.55 counts s _1 (~6. 1 mCrab) in the 20-40 keV band. Timing analysis performed on the ISGRI data 
identifies an orbital period of 4.5696 ± 0.0009 days and gives an ephemeris of mid-eclipse as, T m u = 52931.37 ± 0.04 MJD. 
Re-analysis of archival BeppoSAX data has provided a mass function for the donor star, f(m) = 16 ± IMq and a projected 
semimajor axis of a v sin; = 68 ± 1 lt-s. We conclude that the donor is an OB-supergiant with a mass of 18.8-29.3 Mq and 
a radius of 15.0-23.4 Rq. Spectra obtained by XMM-Newton and ISGRI indicate a high hydrogen column density of N// = 
6.8 x 10 22 cnT 2 , which suggests intrinsic absorption. The source appears to be a high mass X-ray binary with the neutron star 
emitting X-rays through wind-fed accretion while in an eclipsing orbit around an OB-supergiant. 

Key words. X-rays: individuals: IGR J18027-2016, SAX J1802.7-2017 - X-rays: binaries - pulsars: general - Gamma-rays: 
observations 



1. Introduction 

High mass X-ray binaries (HMXBs) were discovered as bright 
X-ray sources in the 1970s. The systems comprise a compact 
object orbiting a massive OB class star. The compact object 
can be either a neutron star or black hole and is a strong X-ray 
emitter via the accretion of matter from the OB companion. 

HMXBs can be split into two individual classes: (/) Be/X-ray 

J II i 

binaries; (z'z) Supergiant X-ray binaries (Corbet, 1986). The ma- 
jority of HMXBs belong to the first class (-80%) (K aperet all 
2004|- In Be systems, the compact object is a neutron star 
and is typically in a wide, moderately eccentric orbit and 
spends little time in close proximity to th e dense circumstel- 
lar disk surrounding t he Be companion, dNegueruela & Coel 
l2002HLiuetail E()00). As a result, they are often X-ray tran- 
sients as the accretion principally occurs when the compact ob- 
ject passes through the Be-star disk. 

In the supergiant systems, the compact object is in orbit 
around an OB-supergiant. Accretion can be powered solely 
by the stellar wind from the supergiant and/or can be through 
Roche-lobe overflow. Accretion from the stellar wind results in 
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an X-ray luminosity of 10 35 — 10 36 erg s . When the star fills its 
Roche lo be a much higher X -ray luminosity is achieved, ~10 38 
erg s -1 ( Kap er et all 1200 4). The orbital periods of the super- 
giant systems are typically shorter (1.4-41.5 days) than the 
Be-star systems (12.7 - 262 days) and the orbits typically more 
circular. 

HMXBs are young systems, and the compact object is of- 
ten a highly magnetized neutron star. Hence, many systems are 
seen to host an X-ray pulsar. The spin period of the pulsar and 
the orbital period of the system correlate well allowing the dif- 
ferentiation between OB-supergiant and Be-star systems and 
between underfilled and filled Roche-lobe supergiant systems 
dCorbeil 19861) . 

Since its launch in 2002, INTEGRAL (the INTErnational 
Gamma-Ray Astrophysics Laboratory) has been performing a 
regular survey of the Galactic plane and a deep exposure of the 
Galac tic Centre as part of its Core Programme JWinkler et all 
2003). In the course of these observations a number of new 
unidentified sources emitting in the soft y-ray region have been 
discovered. The first IBIS/ISGRI soft y-ray galactic plane sur- 
vey catalogue (Bi rd et all |2004) reports the discovery of 28 
objects of unknown classification exhibiting persistent hard 
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X-ray emission. One of the new unidentified sources ob- 
served in the first year of operations of INTEGRAL was IGR 
J18027-2016. This new source has been independently re- 
ported in both core prog ramme GCDE (G alactic Centre Deep 
Expo sure) observations dBird et all |2004) and guest observer 
data dRevnivtsev et alll2004l) ~ 

IGR J18027-2016 is spatially associated with a recently 
identified X-ray pulsar SAX J1802.7-2017 dAu gello et all 
l2003l) . lLutovinov et al.1 ll200.5l) report that the 18-60 keV spec- 
trum of IGR J 18027-20 16 can be described by a power law 
with an exponential cut-off at high energies of E cu , = 18 
keV. We report here the results of the timing analysis of IGR 
J 18027-20 16 using INTEGRAL data, a re-analysis of archival 
BeppoSAX data, and the spectral analysis of XMM-Newton 
data. 



2. Observations 

2.1. BeppoSAX data 

lAugello et al .(2003), hereafter Au03, reported the discovery of 
the serendipitous source SAX J1802. 7-2017, identified as anX- 
ray pulsar. This result was from BeppoSAX archival observa- 
tions of the GX 9+1 field performed from September 16 2001 
(02:01:30.0 UTC) to September 20 2001 (03:00:08.5 UTC). 
The position of SAX J 1802.7-20 17 was given as R.A. (2000.0) 
= 18 /, 02 m 39.9 i and Dec. = -20°17'13.5" with a positional un- 
certainty of 2'. A pulse period of 139.612 s and an orbital pe- 
riod of ~4.6 days was reported. Tablef^shows all of their mea- 
sured system parameters. The total duration of the observation 
~4.0 days, was such that the data did not cover a complete or- 
bit. As a result the measurement of the orbital period through 
the delays in pulse arrival times was assocated with a 9% error 
which contributed to other system parameters having errors of 
~30%. The possibility that this was an eclipsing system was 
also proposed. 



2.2. INTEGRAL data 

IGR J18027-2016 has been observed on multiple occasions 
during the core programme by the IBIS detector ISGR I 
[INTEGRAL Soft Gamma-Ray Imager, dLebrun et alll2003h l. 
IBIS generates images of the sky with ~30° field o f view in 
the energy range 15-1000 keV dUbertini et al.ll2003h . IBIS ob- 
served IGR J18027-2016 over three epochs: 28 Feb. - 28 Apr. 
2003; 27 Sep. - 12 Oct. 2003; 17 Feb. - 20 Apr. 2004. ISGRI 
images were generated for each selected pointing in 10 narrow 
energy bands with the ISDC Offline Scientific Analysis (OSA) 
software version 4. 1 . The individual images were combined to 
produce mosaics of the reg ion in broade r energy bands using 
the system described in Bird et alJ (200l. 

Fig. G] shows the IBIS/ISGRI 20-40 keV image of the re- 
gion around IGR J18027-2016. The atoll source, GX 9+1 is 
also clearly visible -22' from IGR J18027-2016. IGR J18027- 

2016 appears as a 30.5cr detection in the 20-40 keV band at 
R.A. (2000.0) = 18*02 m 46.F a nd Dec. = -20°17 '37.1" with 
a positional uncertainty of ~1 ' faros etallEool . This posi- 
tion is consistent with the reported position of SAX J 1802.7- 

2017 (Au03). The source flux in the 20-40 keV energy band is 
F20-40fev = 0.58 + 0.02 counts sec -1 which corresponds to ~6.4 
mCrab. 



GX9 + 1 



IGR J18027- 




Table 1. Estimated phase delay fit parameters and system pa- 
rameters of SAX J1802. 7-2017 (Au03). [ao, a\, and B are phase 
delay fit parameters defined in Section 3.2.] 



Parameter Value 



ao 
«i 
B 

P orb 

a x sin ( 

T„/2 
* pulse 

KM) 

M c 
R c 



04 +U12 

0.02 ± 0.03 days" 1 
50 +0 07 
4.6+^ days 
70+}° lt-s 

52168.22!° |° MJD 
139 612 +0 ° 06 s 
17 ± 5 M Q 
0.64 ± 0.14 rad 

> 11 M Q 

> 14i? 



Fig. 1. ISGRI 20-40 keV image of the field containing IGR 
Jl 8027-2016 and GX 9+1. The crosses correspond to the 
sou rce positions; the larg er circle is the BeppoSAX error cir- 
cle dAueello et all 120031): the smaller circle is the ISGRI 90% 
error circle dGros et alll2003l 



2.3. XMM data 

IGR J18027-2016 was observed by XMM-Newton for -12 
ksec, from 06:54:40.0 - 09:44:38.0 UTC, 2004 April 06. A sin- 
gle X-ray source was found within the ISGRI err or circle in the 
EPIC PN and MOS cameras dStriider et al 11200 it iTurner et all 
l200ll) . The XMM Science Analysis System (SAS), software 
version 6.0, was used to extract light curves and spectra for IGR 
J 18027-20 16. Standard SAS tools were used to calculate the 
instrumental response and the effective area for the extracted 
spectra. 
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Fig. 2. 20-80 keV light curve of IGR J18027-2016 covering the three observed epochs. The light curve has been rebinned to 2.5 
day averages. 



3. Temporal Analysis 

3.1. ISGRI 

The flux of IGR J18027-2016 was extracted from each ISGRI 
pointing image to generate light curves in different energy 
bands, the rebinned 20-80 keV light curve spanning ~417 days 
is seen in Fig. [2] The source is not especially strong above 
40-50 keV and to minimize systematic effects we searched 
the 20-40 keV light curve for periodicities using the Lomb- 
Scargle periodogram method by means of the fast implementa- 
tion of Press & Rvbicki (1989). The resulting power spectrum 
is shown in Fig. [3] The peak power of 52.58 corresponds to a 
freq uency of 0.2 1 88 days" 1 . The error on the angular frequency 
llHorne & BaliunasU 19861) is: 



5a> 



3ncr N 
l^NTA 



(1) 



where cr N 2 is the variance of the noise, N is the number of 
data points, T is the total length of the data set and A is the 
amplitude of the signal given by: 



A = 2 



N 



(2) 



where zo is the Lomb-Scargle power and cr s 2 is the variance of 
the light curve. 

This frequency corresponds to a period of 4.570 ± 0.003 
days which is consistent with the orbital period indirectly mea- 
sured by Au03. 

The significance of this peak was confirmed by applying a 
randomization test. The time stamps for each flux measurement 
were randomly reordered and the periodogram of the resulting 
light curve generated. We simulated 5000 light curves in this 
fashion. The mean maximum power found in the randomized 
data sets was 10.6 and highest power recorded was a single 
instance of 19.5. It is immediatley obvious that a peak of power 
~53 is extremely unlikely to occur by chance and hence our 
detection is highly significant. 

Upon close examination of the periodogram it can be seen 
that there is some aliasing around the frequency of the detected 
period. This raises the question of whether we have selected the 
correct peak and hence the correct frequency/period. The much 




1.2 0.3 
Frequency (days -1 ) 



Fig. 3. Lomb-Scargle periodogram generated from the 20-40 
keV light curve of IGR J18027-2016. 
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Fig. 4. Top: Zoom on the peak of the periodogram of IGR 
Jl 8027-2016, shown in Fig. [3] Bottom: Lomb-Scargle peri- 
odogram generated from a sine wave with a period of 4.570 
days and sampled identically to IGR J18027-2016. 
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higher power of our selected frequency is sufficient that we can 
be confident that this is the true frequency of the system. To 
investigate the orgin of the aliasing pattern a pure sine wave 
of our chosen frequency was sampled in an identical fashion 
to that of the original data set. The corresponding periodogram 
(Fig. |3 can be seen to exhibit broadly the same aliasing pattern 
as observed in the data. The peak of maximum power corre- 
sponds to the frequency of the input sine wave and the aliases 
are spaced in an identical fashion to that of the original data 
set. Also, as discussed in section 3.2, folding the ISGRI light 
curve on a period of 4.570 days indicates an eclipse of the sys- 
tem exactly when we would expect it to be observed. Hence, 
we are confident that the peak of maximum power in our data 
corresponds to the true frequency of the system. 

3.2. BeppoSAX 

We used the orbital period as measured by ISGRI to re-analyse 
the BeppoSAX data presented by Au03. Originally, phase de- 
lays in the pulse period of the pulsar measured by SAX were 
fitted without any prior knowledge of the orbital period. We 
performed an identical analysis of the phase delays assuming 
an orbital period of 4.570 days. 

The modulation of the phase delays can be explained by 
the propagation delays resulting from the orbital motion of the 
neutron star around its companion star. The phases were fit with 



A(p = qq + a 1 1„ + B cos 



2n(t„ - T„n) 



orb 



(3) 



where t„ is the arrival time of the nth pulse, T„/2 is the 
epoch of NS superior conjunction of the NS, the linear term 
ai = AP pulse /P 2 pulse and P orh = 4.570 + 0.003 days. The^ of 
the fit was 1.3, the fit parameters are listed in Table 13 and the 
fit is shown in Fig.|5] As a x sm.i = P pl ,i seB-> then £z A -sinz ~ 68 ± 1 
lt-s; Ppuhe ~ 139.612 s (Au03). All of the fit parameters are con- 
sistent with those originally measured by Au03, however they 
are considerably better constrained due to the independent and 
precise measurement of P orc , by IBIS/ISGRI. 



Table 2. Phase delay fit parameters 



Parameter 


Value 


ao 


0.03 ± 0.07 




0.02 ± 0.03 days" 1 


B 


0.49 ± 0.003 


T n /2 


52168.26 ± 0.04 MJD 



Using the measured orbital period and the epoch of superior 
conjunction of the NS, T„p_, obtained in the fit, the IBIS/ISGRI 
light curves were phase-folded in the 20-30 keV and 30-40 keV 
energy bands. These are shown in Fig. [5] and clearly show a 
sharp eclipse as suggested by Au03. However, the eclipse ev- 
ident in Fig. [6] appears assymmetic and indicates that we may 
be partially resolving the eclipse suggesting that the hard X-ray 
emission may not exclusively originate from a point source. 



f-rti W 



TIME(doys) - 52168.1 1 MJD 



Fig. 5. Phase delays as a function of time as measured by 
BeppoSAX. The dashed line is the best-fit function (Eq.Q). 



o.o 

0.60 
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20-30 keV 
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Fig. 6. Folded light curves in the 20-30 and 30-40 keV energy 
bands of IGR J18027-2016. The folding is performed for an 
orbital period of 4.570 days, and the zero epoch was assumed 
at the superior conjunction of the NS,T ff /2 = 52168.26 MJD . 



The fuller phase coverage of the IBIS/ISGRI data unambigu- 
ously confirms the orbital period and the presence of an eclipse 
in the system. The time of mid-eclipse occurs at the time of 
superior conjunction of the NS, T„/2. 
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3.3. Refining the orbital period & an accurate 
ephemeris 

Taking the zero epoch as the time of superior conjunction of 
the NS, T„ /2 = 52168.26 ± 0.04 MJD and the orbital period, 
Porb - 4.570 ± 0.003 days the number of orbital cycles which 
have occured by the time of the last ISGRI observation is E 
= 207.2 ±0.1. Hence, over this timescale, we can confidently 
identify which cycle we observe. This is borne out by Fig. [6] 
as the eclipse is evident very near zero phase when the ISGRI 
data is folded with the BeppoSAX zero epoch. From the folded 
ISGRI light curves the time of an eclipse around the mid point 
of the ISGRI observation is calculated, T eclipse = 52931.4 ± 0.2 
MJD. This would correspond to the beginning of the 167 th or- 
bital cycle. Linearly fitting between these two times of eclipse 
measurements yields a refined orbital period estimate of P or b = 
4.5696 ± 0.0009 days. This allows the calculation of an accu- 
rate mid-eclipse ephemeris from: 



Tmid = Tq + P orb x E 



(4) 



where To is the time of the 1" eclipse seen by BeppoSAX 
and E is the cycle number. This gives an ephemeris of, r m ,y 
= 52931.37 ±0.04 MJD. 

3.4. XMM 

The timing mode of the PN camera all ows the collection o f tim- 
ing data with a resolution of 30 fis ( Striid er et all 12001 ). The 
arrival time of all events were corrected to the solar system 
barycentre, using the SAS task "barycen". A 2-10 keV band 
light curve consisting of 1 second bins was generated for the 
duration of the observation. The Lomb-Scargle periodogram 
method performed on the ISGRI data was applied to search for 
periodicites. The resulting power spectrum is shown in Fig.Q 
A peak of power ~504 is clear at a frequency of 7.17 mHz, 
the second harmonic of this frequency is also clearly visible 
with a power of ~220. This corresponds to a pulse period of, 
P pulse = 139.47 ± 0.04 seconds. This is consistent with the 
range of pulse periods (139.44 - 139.86 s) detected by Au03. 




0.008 0,010 
Frequency (Hz) 




Spin Phase 



Fig. 8. Pulse period phase folded light curve in the 2-10 keV 
energy band. The folding is obtained for a pulse period of 
139.47 seconds and the zero epcoh is taken at the superior con- 
junction of the NS, T„ i2 = 52168.26 MJD. 



At the time of the XMM observation we expect the system to 
be at a phase of ~0.2 in its orbital cycle. Hence, the pulse period 
should be doppler shifted by the orbital motion. Correcting the 
measured pulse period to the rest frame of the pulsar yields, 



pulse 



= 139.61 ± 0.04 seconds. The XMM measurement is 



Fig. 7. The Lomb-Scargle periodogram generated from the 
XMM EPIC 2-10 keV light curve of IGR J18027-2016. 



consistent with the pulse period, P pu i se = 139.612 ± 0.006 sec- 
onds, measured by BeppoSAX. 

Folding the XMM data on this pulse period yields the pulse 
period phase folded light curve seen in Fig. [8] The main pulse 
occurs around phase 0.63. A secondary peak is visible around 
phase 0.05. The pulse fraction, (I„ mx - 7„„„)// mflA , where I max 
and I m j„ are the maximum and minimum count rates, is ~68% 
± 10%. The pulse period phase folded light curve matches well 
with the folded light curves of Au03, the pulse fractions are 
consistent and the locations of the primary and secondary peaks 
also agree. 

4. Spectral Analysis 

The ISGRI spectrum was extracted from narrow energy band 
mosaics and hence represent the long term average spectrum 
of the source. Phase resolved spectra were extracted from the 
XMM-EPIC data set. Using the pulse profile (Fig. |SJ a spec- 
trum of the source during the primary pulse and a spectrum 
outside of the primary pulse were generated. Au03 reported that 
the low statistics of the BeppoSAX data did not allow for an ac- 
curate spectral analysis. Additionally, the BeppoSAX PDS data 
is difficult to analyse as the data is contaminated by the bright 
atoll source, GX 9+1 which also appears in the field of view. 

The EPIC PN and ISGRI average spectra were simultane- 
ously fit, in XSPEC vl 1.3, by a photoelectric absorbed, broken 
power law model with a single Gaussian emission line (gauss + 
zvphabs*bknpower). A broad Gaussian line is used to account 
for soft X-ray residuals that mostly concentrate at low energy. 
This soft excess contributes 7 xl0~ 14 ergs cirT 2 s _1 below 3 
keV. Elemental abundances of the spectral model, excluding 
Fe, are frozen to solar quantities and the redshift is frozen to 0, 
all other parameters are free. The fit has;^ = 1.29. The best fit 
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channel energy (keV) 



Fig. 9. The phase resolved XMM (upper left: EPIC PN spectrum of primary pulse. Lower left: EPIC PN spectrum outside of the 
primary pulse) and INTEGRAL photon spectrum (right: ISGRI spectrum)of IGR J18027-2016 with the best fit model. 



parameters are listed in Table[3] The broken power law indices 
were found to be Ti = 0.79 ± 0.08 and T 2 = 3.1 ± 0.1, with a 
break at Egreak = H-3 ± 0.6 keV. A column density of Nh = 
6.8 ±1.0 xlO 22 crrT 2 is found compared to the expected galac- 
tic column density of 1.0 xlO 22 cirT 2 . The Fe edge is ~3 times 
that of solar abundances, this may indicate that large amounts 
of Fe are present or that the Nh is underestimated. The shape 
of the spectrum is typical of X-ray pulsars which are character- 
ized by a flat power law with a slope of between and 1 up to a 
high energy cu t-off at 10-20 keV a fter which the slope becomes 
much steeper ( Wh ite et all 1 1 995 ) . Fitting Fe lines at 6.45 and 
7.1 keV improves the overall fit, however the F test indicates 
that they are significant at only the 95% confidence level. The 
ISGRI spectrum flux is, F2o-iooteV = 5.9 xl0~ n ergs crrT 2 s , 
the XMM unabsorbed fluxes are, Fa-iofev = 8.9 xl0~ n ergs 
cm' 2 s _1 within the primary pulse, and ¥i-\okeV - 5.1 xl0~ u 
ergs cirT 2 s _I outside the primary pulse. 

The shapes of the phase resolved XMM spectra are consis- 
tent with each other. The hydrogen column density and contin- 
uum parameters of the model, excluding the normalization, do 
not vary with pulse phase indicating that the absorption is not 
associated with the neutron star accretion column. 

It is possible to fit the data with a Comptonization model or 
a cut-off power law, however to do so requires the introduction 
of a particularly small normalization constant to the spectral 
model, C/sgr/~0.2, to account for cross-calibration uncertain- 
ties of the ISGRI and EPIC observations. This would indicate 
that XMM had observed IGR J18027-2016 while it was in a 



Table 3. Best fit spectral parameters (90% confidence) of the 
spectral model discussed in the text and simultaneously fit to 
the XMM-EPIC and INTEGRAL-ISGRI spectra. 

Parameter Value 

X 2 /v 164.57/128 

N H 6.8 + 1.0 x 10 22 cirT 2 

r, 0.79 ± 0.08 

Eateak 11.3 ± 0.6 keV 

T 2 3.1 ±0.1 

Bkn power law I lkeV 4.9 ± 1.0 x 10 -3 ph keV~' cnT 2 s" 1 

Fe abundance 3.1 ± 0.8 (of solar) 



high state when compared to the average ISGRI flux. However, 
the ISGRI light curve (Fig.0 at the time of the XMM observa- 
tion show no indication of a large change in flux of the source. 
The fitting of system Nh is highly dependent upon the spectral 
shape of the continnum assumed in the spectral model, a softer 
spectrum (e.g. the comptt model describing Comptonization of 
soft photons in a hot plasma) results in an Nh twice as large as 
listed in Table|3] 

5. Discussion 

We have directly and precisely measured the orbital period 
for IGR J18027-2016 of 4.5696 ± 0.0009 days. This is con- 
sistent with the implied orbital period of SAX J1802.7-2017, 
and as it is the only known X-ray source within 10' accord- 



A. B. Hill et al.: The 1-50 keV spectral and timing analysis of IGR J18027-2016 



7 




2\) 25 3D 35 HI 

iZjperqant moss (Soar Messes) 



Fig. 10. The derived mass-radius relationship of the supergiant (SG) in IGR J 18027-20 16 (assuming a NS mass of 1.4 Mq) is 
plotted against the SG roche lobe radius and the NS-SG separation. The lighter shaded area corresponds to the region in which 
the SG radius is smaller than its roche lobe radius. Also plotted are the masses and radii of a number of known OB supergiants. 



ing to the SIMBAD/NED database the two sources are con- 
firmed as the same source. Additionally, we have observed that 
this is an eclipsing system confirming the hypothesis of Au03. 
The ephemeris of the mid-eclipse time is measured as, T m M = 
52931.37 ± 0.04 MJD. The -139 s pulse period measured by 
BeppoSAX and XMM combined with the 4.5696 day orbital 
period, situates I GR J1 8027-2016 in the region of the Corbet 
diagram JCorbetl [l986) occupied by underfilling Roche lobe 
supergiant systems . Hence, the system appears to be a neutron 
star that is accreting material through the stellar wind of the 
companion supergiant star. 

The ISGRI 20-80 keV light curve, shown in Fig. indi- 
cates that this is a weak, persistent s ource with no evidence 
of flaring in any of our observations. iLutovinov et al] (|2005) 
also reported that the 18-60 keV fluxes observed by IBIS in 
September 2003 and March 2004 of IGR J18027-2016 were 
nearly constant and equal. The absence of strong hard X- 
ray variability contrasts with what has been observed in other 
HMXBs in which flares by a factor of 10 have been observed. 
This may indicate that the stellar wind is rather homegenous at 
least on time scales longer than a few days. The persistent na- 
ture of the source suggests that the NS is continuously accret- 
ing, and, combined with the sinusoidal modulation of the pulse 
arrival times this implies that the NS is in an approximately 
circular orbit about its companion. This is consistent with the 
eccentricity upper limit of e < 0.2 given by Au03. The absence 



of apparent transient behaviour and circular orbit further sug- 
gests that this is a supergiant system and not a Be X-ray binary 
Liu et allEoOOl) . However, the fact that IGR J18027-2016 has 
not been reported by previous X-ray telescopes may indicate 
some level of transient behaviour. Notably there was no de- 
tection by the ROSAT mission. Transient behaviour would in- 
dicate the absence of Roche lobe overflow but would suggest 
that this was an atypical Be X-ray binary as its orbital and pulse 
period do not locate it in the expected region of the Corbet di- 
agram and it has a very low eccentricity orbit. 

During the XMM observation the X-ray flux increased by 
50% with constant pulse fraction. Spectra derived for both peri- 
ods do not show any significant variation in the spectral param- 
eters. This is notable as one could expect to see a significant in- 
crease in Nh for large flux differences. The observed properties 
of IGR J 18027-20 16 appear to be very stable on timeescales of 
a few days, as seen by ISGRI, and 12 ksec, as seen by XMM. 

The mass function of the system donor is: 

Mr sin 2 / / 2n \ 2 (a x sin/) 3 

where M c is the supergiant companion star mass, q — 
is the ratio of the NS mass to the supergiant companion mass 
and G is the gravitational constant. 

Fig. |6] clearly shows an eclipse occuring in the 20-30 and 
30-40 keV energy bands. Future observations of the ingress 
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and egress of the NS eclipse may provide information regard- 
ing the donor star's atmosphere and indicate any evidence of 
extended high energy emission in the region around the NS. 
There are two methods to determine the half angle of eclipse, 
9 e , which we define as the half width half height of the eclipse. 
Firstly, if we assume a cicular orbit, the measured duration of 
the eclipse corresponds to a half angle of eclipse of e - 0.66 ± 
0.2 radians in both energy bands. Additionally, the mid-eclipse 
time corresponds to the epoch of NS superior conjunction of 
the NS, T„/2, derived from our fit. This corresponds to the pe- 
riod in the BeppoSAX data when the count rate was consistent 
with the background count rate. Taking the time at which the 
eclipse ends as 52,168.69 MJD (Au03), then the half-angle of 
the eclipse corresponds to 9 e = 0.591 + 0.09 radians. Taking 
a weighted mean of the half angle of eclipse gives 6 e = 0.61 
± 0.08 radians. Assuming we have a point source and spher- 
ical geometry the duration of the eclipse is related to Rc, the 
companion star radius, th e inclination of the system and a, the 
separation of the stars by (Nagase, 1989): 



Table 4. System parameters of IGR J18027-2016 



— = (cos 2 i + sin 2 sin 2 G e ) 
a 



1/2 



(6) 



Assuming that the NS has the canonical mass of 1.4 Mq 
then for any mass of the companion star the inclination of the 
orbit can be dervied from the mass function (eq.|5j. Assuming 
that we have a perfectly circular orbit then a x ^a, the phase de- 
lay modulation yields a x sin i and hence from eq.|5]there is an 
empirical relationship between the companion star mass and ra- 
dius. An analytic approximation for the companion star Roche 
lobe is: 



R 



Roche 



a 0.38 + 0.2 log 



M c 

M, 



A'.V 



(7) 



■ 1 1 i 

(Bowers & Deeming, 1984) (accurate < 1%) 

From its low luminosity and its position in the Corbet dia- 
gram, we expect the supergiant companion star to be underfill- 
ing its Roche lobe hence the Roche lobe radius gives an upper 
limit to the allowed companion star radius, mass and system 
inclination. The empirical Mc-Rc relation and the associated 
RRoche are plotted in Fig.[K)]along with the neutron star - super- 
giant separation. The lighter shaded region towards the lower 
left of Fig. 1101 indicates the allowed system parameters such 
that our measurements represent an underfilling Roche lobe su- 
pergiant. Also plotted a re the masses and radii of a number of 
know n OB supergiants dWilson & DopitaLll985tlHerrero et all 
Il992h . 

From Fig.EI|we can see that the properties of the compan- 
ion star are: Mc ~ 21Mq; Rc ~ 19Rq. The complete range of 
these values is given in Table|4]with the other parameters of the 
system. These mass and radius values are consistent with late 
09 - early Bl supergiant stars. 

Measurement of changes in the NS pulse period is an inde- 
pendent indicator of the accretion mechanism operating in the 
system. If the NS is accreting matter from the donor star via 
Roche lobe overflow the NS should show evidence of spinning 
up. The BeppoSAX and XMM measurements of the spin pe- 
riod as 139.612 ± 0.006 s and 139.61 + 0.04 s respectively, are 



Parameter 



Value 



Porb 

P, mhe (BeppoSAX) 
P pulse (XMM) 
a x sin i 

T„/2 
Pmid 

f(M) 

e e 

M c 
Rc 
a 
i 



4.5696 ± 0.0009 days 
139.612 ± 0.006 s 
139.61 ± 0.04 s 
68 ± 1 lt-s 

52168.26 ± 0.04 MJD 
52931.37 ±0.04 MJD 
16 ± 1 M Q 
0.61 ± 0.08 rad 
18.8 < M < 29.3 M Q 
14.7 < R < 23 A Rq 
28.4 < a < 35.8 Rq 
59° < i < 88° 



Table 5. Theoretical P pu i se and distance for different source lu- 
minosity 



Luminosity, (ergs cm 2 s 1 ) 


10 35 1 36 1 37 1 38 


Ppuise, (s yr" 1 ) 
Distance, (kpc) 


0.02 0.2 1.4 9.7 
3.1 9.7 31 421 



separated by ~2.5 years and indicate no evidence of a change 
in spin period. Taking the lcr limits of these measurements we 
find an upper limit to the spin-up rate of, P nll i^~0 . 02 s y r -1 . 
The theoretical spin-up rate is given by (Fra nk et all ll992): 



^^8x 10- 5 M- 3/ X /? Ll^lo^Pp^yr- 1 



pulse 



(8) 



where Mi is the NS mass in Mo, Re is the NS radius in cm, 
L37 is the NS accretion luminosity in units 10 37 ergs cirT 2 s _1 , 
^30 is the magnetic moment in units 10 30 G cm 3 and I45 is the 
moment of inertia (~MR 2 ) of the NS in units 10 45 g cm 2 . For a 
typical neutron star with a magnetic field strength, B~10 12 G, 
A*30 ~ I45 ~ R6 ~ 1 ■ The theoretical spin-up of a neutron star ac- 
creting via Roche lobe overflow with a luminosity in the range 
10 35 -10 38 ergs cirT 2 s' 1 is shown in Table|5]with the associated 
source distance for each luminosity. If Roche lobe overflow is 
the source of accretion then a luminosity of 10 38 e rgs cm' 2 s -1 , 
close to the Eddington luminosity, is expected (K aoer et all 
2004). Table |5]illustrates that this would locate the source out- 
side of the galaxy and would result in a clear spin-up of the 
pulsar of ~24 seconds between the BeppoSAX and XMM ob- 
servations, which is not seen. The limit of P p „i se ~0.02 s yr' 1 , 
would imply that the source has a luminosity of 10 35 ergs cm -2 
s' 1 which is very much lower than expected from a Roche lobe 
overflow system but would be typical of a wind accretion sys- 
tem. Hence, this is another indication that the donor star has 
not filled its Roche lobe and that the pulsar is powered en- 
tirely by wind fed accretion and thus lending further credence 
to the constraints on the donor star given by the shaded region 
of Figure [TO] 

The soft excess (< 3keV) evident in the XMM-EPIC data 
is identical in the on and off pulse spectra implying that the 
soft excess may not be pulsed. Performing timing analysis on 
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the events below 3 keV does not show any period, however, the 
event rate is not high enough to reach a definitive conclusion. 
Soft X-ray excess have been found in other HMXB sources, 
e.g. Vela X-l, and have been explained by scattering or partial 
ionisation of the stellar wind. 

The equivalent width of the 6.4 keV Fe line has a 3 sigma 
upper limit of 40 eV and 25 eV for the off pulse and on pulse 
spectra respectively. Such values are consistent with Nh < 5.0 
x 10 22 cirT 2 for a spherical geometry. The non-detection of 
the Fe line is therefore consistent with the absorbing column 
density which is oberserved. 

The spectral fit of the EPIC and ISGRI spectra indicate a 
column density of Nh = 6.8 ± 1.0 x 10 22 cm~ 2 . This is much 
greater than the expected line of sight column density, Nh = 
1.0 x 10 22 crrT 2 , and suggests that the absorption is intrinsic 
to the source and not the line of sight. This could be explained 
by the stellar wind from the OB supergiant accreting onto the 
neutron star forming a dense spherical shell. The low statistics 
of the BeppoSAX data used by Au03 did not allow an accurate 
spectral analysis. They inferred a luminosity of ~5.6 x 10 35 
ergs s _1 over the 1.8-10 keV range by fitting their data to a 
power law with F = -0.1 ± 0.1 and N H = 1.7^ x 10 22 cirT 2 , 
and assumed a source distance of 10 kpc. The XMM spectral 
fit indicates that the column density is much higher and hence 
the source luminosity was previously underestimated. 

Assuming that IGR J 18027-20 16 is 10 kpc away the 2-10 
keV flux during the primary pulse of F2-\ok e v = 8.9 xlO 11 ergs 
cirT 2 s _1 gives an estimate of the source luminosity of ~1.1 
x 10 36 ergs s -1 , a typical luminosity for a non-Roche lobe filling 
Supergiant HMXB ( Ne gueruelal 12004). A useful comparison is 
provided by 4U 1538-52, which is a typical HMXB accreting 
from a stellar wind and shows similarities to IGR J18027-2016. 
4U 1538-52 has a long spin period of ~529 s, an orbital pe- 
riod of 3.75 days and undergoes a well de fined X-ray eclipse . 
The companion is a 19.9Mq B0 star jRevnolds etafl 1 19921) 
and the X-ray flux has been estimated as ~ 2 x 10 36 ergs s _1 for 
a distance of 6.4 kpc (Robb a et all 1200 II) . Hence the luminos- 
ity of IGR J18027-2016 appears to be typical of wind accreting 
HMXBs assuming it has a distance which approximates to the 
far side of the Galactic Centre. 

6. Summary 

We presented an analysis of INTEGRAL, BeppoSAX and 
XMM-Newton observations of IGR J18027-2016. We con- 
clude that the source is an eclipsing HMXB system comprising 
an X-ray pulsar accreting matter from the stellar wind of a late 
O - early B supergiant. The source is persistent and has a high 
intrinsic photoelectric absorption and is similar in this nature 
to a n umber of other new INTEGRAL sources ( I Walter et all 
120031) . IGR J18027-2016 has the properties of a typical wind 
accreting HMXB system. 
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